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Abstract 

In this paper, we propose a resource allocation algorithm for the downlink of sectorized two- 
dimensional (2-D) OFDMA cellular networks assuming statistical Channel State Information (CSI) and 
fractional frequency reuse. The proposed algorithm can be implemented in a distributed fashion without 
the need to any central controlling units. Its performance is analyzed assuming fast fading Rayleigh 
channels and Gaussian distributed multicell interference. We show that the transmit power of this simple 
algorithm tends, as the number of users grows to infinity, to the same limit as the minimal power required 
to satisfy all users' rate requirements i.e., the proposed resource allocation algorithm is asymptotically 
optimal. As a byproduct of this asymptotic analysis, we characterize a relevant value of the reuse factor 
that only depends on an average state of the network. 

I. Introduction 

We address in this work the problem of resource allocation (power control and subcarrier assignment) 
for the downlink of sectorized OFDMA networks impaired with multicell interference. A considerable 
research interest has been lately dedicated to this problem since the adoption of OFDMA in a number of 
current and future wireless standards such as WiMax and 3GPP-LTE. In principle, the problem of resource 
allocation should be jointly solved in all the cells of the system. In most of the practical situations, this 
optimization problem is difficult to solve. Therefore, most of the related works in the literature focus on 
the single cell case {e.g., [l]-[6]). Fewer works address the more involved multicell allocation problem. 
In this context, we cite [7]-[ll] in the case of perfect CSI at the transmitters side, and [12], [13] in the 
case of imperfect CSI. In [12], [13], all the available subcarriers are likely to be used by different base 
stations and are thus subject to multicell interference. In such a configuration, interference may reach 
excessive levels, especially for users located at cells borders. 



September 14, 2010 



DRAFT 



2 



Similarly to [12], [13], we assume in this paper that users' channels undergo fast fading and that the 
CSI at the base stations is limited to some channel statistics. However, contrary to these two works, we 
consider that a certain subset of subcarriers is shared orthogonally between the adjacent base stations 
(and is thus "protected" from multicell interference) while the remaining subcarriers are "non protected" 
since they are reused by different base stations. This so-called fractional frequency reuse (or FFR) is 
recommended in a number of standards e.g., in [14] for IEEE 802.16 (WiMax) [15], as a way to avoid 
severe inter-cell interference. The ratio between the number of non protected subcarriers and the total 
number of subcarriers is generally referred to as the reuse factor and is denoted in the sequel by a. 

Few works in the literature (we cite [16], [17], [18] without being exclusive) have addressed the 
problem of resource allocation for FFR-based OFDMA networks, and none of them fits into the above 
framework which is considered in this paper. The particular problem considered in [16] consists in 
maximizing a system-wide utility function under a power constraint. In this context, the authors propose 
a distributed iterative allocation algorithm that is based on estimating the level of multicell interference 
rather than computing it. Of course, resource allocation schemes that do not resort to such simplifications 
are highly preferable. In the same context, authors of [17] consider the problem of minimizing the total 
transmit power needed to satisfy all users' rate requirements. For that sake, they propose a heuristic 
allocation algorithm without any assessment of its deviation from the optimal solution to the latter 
problem. Moreover, the selection of a relevant reuse factor is not addressed. Finally, authors of [18] 
assume that subcarrier assignment is done separately and in advance. The major drawback of this work 
is thus that joint power control and subcarrier assignment is not addressed. 

In our work, we investigate the problem of power control and subcarrier assignment for the downlink 
of FFR-based OFDMA systems allowing to satisfy all users' rate requirements while spending the least 
possible power at the transmitters' side. In our previous work [19], [20], the solution to this problem is 
characterized in the special case of one-dimensional (1-D) cellular networks where all users and base 
stations are located on a line. Unfortunately, it is much more difficult to characterize this solution in the 
case of 2-D networks. In the present work, our aim is to propose a suboptimal resource allocation strategy 
for these 2-D networks and to study its performance with respect to the above optimization problem. Our 
allocation algorithm assumes that users of each cell are divided prior to resource allocation into two groups 
separated by a fixed curve. The first group is composed of closer users to the base station. These users 
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are constrained to modulate non protected subcarriers and are thus subject to multicell interference. The 
second group comprises the farthest users who are constrained to modulate interference-free subcarriers. 
In order to relevantly select the aforementioned separating curves, we study the limit of an optimal 
solution to the resource allocation problem as the number of users grows to infinity. We show that if 
the curves are set using the results of this asymptotic analysis, then the limit of the transmit power of 
the proposed suboptimal algorithm is equal to the limit of the transmit power of the optimal resource 
allocation. As a byproduct, we are able to determine a relevant value of the reuse factor. Indeed, the 
asymptotic transmit power depends on the average rate requirement and on the density of users in each 
cell. It also depends on the value a of the frequency reuse factor. We can therefore define the optimal 
reuse factor as the value of a which minimizes this asymptotic power. The main contributions of this 
work are thus the following. 

1) A practical resource allocation algorithm that can be implemented in a distributed manner is 
proposed for the downlink of a sectorized OFDMA network assuming fractional frequency reuse 
and statistical CSI. The transmit power of this simple algorithm tends, as the number of users grows 
to infinity, to the same limit as the minimal power required to satisfy all users' rate requirements. 

2) As a byproduct of our study of the above algorithm, we prove that the simple scheme consisting in 
separating users of each cell beforehand into interference-free users (constrained to modulate only 
non reusable subcarriers) and interference users (constrained to modulate only reusable subcarriers) 
is asymptotically optimal. This scheme is frequently used in cellular systems, but it has never been 
proved optimal in any sense to the best of our knowledge. 

3) Finally, a method is proposed to select a relevant value of the reuse factor. The determination of 
this factor is of great importance for the dimensioning of wireless networks. 

The rest of this paper is organized as follows. The system model is introduced in Section II, followed by 
a description of the multicell resource allocation problem in Section III. The proposed resource allocation 
algorithm is presented in Section IV. The relevant choice of the curves associated with this algorithm 
and which separate the two groups of users in each cell is addressed in Section V. Next, the relevant 
selection of the reuse factor is addressed in Subsection V-D. Finally, The relevancy of the proposed 
resource allocation and of our selection of the reuse factor are sustained by simulations in Section VI. 
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II. System Model 



Consider the downlink of a sectorized OFDMA cellular network composed of hexagonal cells. Each 
cell in the system is divided into three 120° sectors. In this paper, we restrict ourselves to the case of 
three interfering sectors of three adjacent cells, say cells A,B,C (see Figure 1). In the more general 
case of networks with more than three cells, our results hold provided that the interference generated by 
farther base stations can be neglected. Generally, this assumption is only valid as a first approximation. 
However, it allows for an essential reduction of the dimensionality of the multicell resource allocation 
problem. In the sequel, we assume that the considered sectors of cells A,B,C have the same surface 
and we denote by K"^, , K'^ their respective number of users. Let K = + + be the total 
number of users and N the total number of available subcarriers. The signal received by user k in cell c 
(c G {A, C}) at subcarrier n G {0 . . . — 1} during the mth OFDM block is given by 



where sj^{n,m) represents the data symbol destined to user k, and where Wk{n,m) is a random pro- 
cess that encompasses both the thermal noise and the possible multicell interference. Random variable 
H^{n, m) stands for the frequency-domain channel coefficient associated with user k in cell c at the 
nth subcarrier and the mth OFDM block. The realizations of this random variable are assumed to be 
known only at the receiver side and unknown at the base station. Random variables {H^{n,m)}n,m are 
Rayleigh distributed with variance pk = E[|i/^(m, which is assumed to be constant w.r.t n and m. 
This holds for example in the case of uncorrelated time-domain channel coefficients. Furthermore, for 
each n G {0 . . . — 1}, random process {Hf, {n, m)}m is assumed to be ergodic. Finally, variance pk 
is assumed to be known at the transmitter side and vanishes with the distance between base station c 
and user k following a given path loss model. We assume that fractional frequency reuse is applied. 
According to this scheme (see Figure 1), a certain subset of subcarriers J C {0 ... A^ — 1} is reused in 
the three cells. If user k modulates a subcarrier n G J, the noise Wk{n,m) includes both thermal noise 
and multicell interference. The reuse factor a is the ratio between the number of reused subcarriers and 
the total number of subcarriers: 



yk{n, m) = H^{n, m)sk{n, m) + Wk{n, m) 



(1) 



a = 



card(J) 
N 
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The remaining (1 — a)N subcarriers are shared by the three sectors in an orthogonal way, such that each 
base station c {c = A, B, C) has at its disposal a subset of cardinality ^-^N . If user k modulates a 
subcarrier n G process Wk{n,ni) will contain only thermal noise with variance a^. Finally, JU CPyiU 
!Pb U 7c = {0, 1, . . . , — 1}. Denote by the subset of subcarriers assigned to user k. We assume that 
Nfc may contain subcarriers from both the "interference" subset J and the "protected" subset ^c- Denote 
by 7^ -^N (resp. 7^ the number of subcarriers assigned to user A; in J (resp. Tc)- In other words, 

7^,1 = card(J n Nfe)/iV = card(y, n Xfc)/iV . 

Parameters 7^ ^ and 7^ ^ generally referred to as sharing factors. We assume from now on that they 
can take on any value in the interval [0, 1] (not necessarily integer multiples of 1/A^). 

Remark 1. Even when the sharing factors are not integer multiples ofl/N, it is still possible to practically 
achieve the exact values of 7^ 7^ 2 ^ifnply exploiting the time dimension. Indeed, the number of 
subcarriers assigned to user k can be chosen to vary from one OFDM symbol to another in such a way 
that the average number of subcarriers in subsets J and 7c is equal to 7^ i-/V and 'j'j^ respectively. 
Thus the fact that 7^ ^, 7^ 2 '^^^ strictly integer multiples of 1 /N is not restrictive, provided that 
the system is able to grasp the benefits of the time dimension. The particular case where the number of 
subcarriers is restricted to be the same in each OFDM block is addressed in Section VI. 

Note that by definition 

7m < a , 7fc,2 < — ^ • 

k=l k=l 

For the sake of readability and compactness of the paper, the above two inequality constraints will be 
written from now on as equalities i.e., we force the whole set of available subcarriers to be fully occupied 
by setting X^aI^i 7fc 1 = a and X]fcLi7fc2 — Indeed, keeping the above constraints as inequalities 
would make the presentation of the final results as well as of the proofs very tedious. 

Recall that in our model, for each user k in any cell c, all channel coefficients H^{n, m) are identically 
distributed on all the subcarriers assigned to this user (the variance pk = E[|//^(?n, is assumed to be 
constant w.r.t n). It is thus reasonable to assume that the base station modulates the subcarriers of each 
user in each one of the two subsets (J and 7c) with the same transmit power. Define P^^ (resp. P^2) 
the power transmitted on the subcarriers assigned to user A; in J (resp. in 7c) i.e., P^i = E[\sk{n,m)\'^] 
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if n S J, i-*^2 = E[\sk{n, m)p] if n G Tc- Parameters {7^ -P^j}i=i,2 will be designated in the sequel as 
the resource allocation parameters. We now describe the adopted model for the multicell interference. 
Consider one of the non protected subcarriers n assigned to user k of cell A in subset J. Denote by a'j. 
the variance of the additive noise process Wk{n, m). This variance is assumed to be constant w.r.t both n 
and m. It depends only on the position of user k and the average powers = Yl,k=i ^ki^ki ^'^^ 
1 ^ki^ki transmitted respectively by base stations B and C on the subcarriers of J. This 
assumption is valid for instance in OFDM A systems that utilize random subcarrier assignment [21]. 
According to this subcarrier assignment scheme, each user k is assigned a subset that is composed by 
randomly selecting card (]\ffc) subcarriers out of the total N available subcarriers. Finally, let cr^ designate 
the variance of the thermal noise. Putting all pieces together: 



¥.\\wk{n,m)\''] 



if n G 

(2) 

I '^1 = + T.c=B,c^[\Hl{nM?]Q\ if^GJ 



where H^{n,m) (c = i?,C) represents the channel between base station c and user k in cell c on 
subcarrier n and OFDM block m. Of course, the average channel gain E [|i7^(n, m)p] depends on the 
position of user k via the path loss model. For instance, if two users k and / of cell A are located on 
the same line perpendicular to the axis BC such that k is closer to base station A, then < af. 

III. Joint Resource Allocation Problem 

Assume that each user k has a rate requirement of Rk nats/s/Hz. Consider the problem of determination 
of the resource allocation parameters for the three interfering sectors. These parameters must be selected 
such that the target rate of each user is satisfied and such that the power spent by the three base stations is 
minimized. Due to the ergodicity of the process {H^{n,m)}m for each subcarrier n, the rate Rj. can be 
satisfied provided that it is smaller than the ergodic capacity Ck associated with user k. Unfortunately, the 
exact expression of Ck is difficult to obtain due to the fact that the noise-plus-interference {wk{n, m)}„ m 
is not a Gaussian process in general. Nonetheless, if we endow the input symbols Sk{n, m) with Gaussian 
distribution, the mutual information between Sk{n, m) and the received signal ?/fc(n, m) in (1) is minimal 
when Wk{n, m) is Gaussian distributed. Therefore, we approximate in the sequel the multicell interference 
by a Gaussian process as this approximation provides a lower bound on the mutual information. Focus on 
cell A and denote by gk,i{Qi -.Qi ), gk,2 the channel Gain-to-Interference-plus-Noise Ratio (GINR) and 
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Gain-to-Noise Ratio (GNR) associated with user k on the subcarriers of subset 3 and "Pa respectively: 

//O-B /-)C\ Pk Pk 
9k,l[Ql , Vl j = — 2 ' Sk,2 = — n • 

The ergodic capacity associated with user k in cell A is equal to the sum of the ergodic capacities 
corresponding to both subsets 3 and Ta- For instance, the part of the capacity corresponding to the 
protected subset Pa is equal to 7^2^ log ^1 + P^2 ^^'' ''"2'"^''^ ^ , where factor 7^2 traduces the fact 
that the capacity increases with the number of subcarriers which are modulated by user k. In the latter 
expression, the expectation is calculated w.r.t random variable \MjJ^i[})\-_ Now, ^^'"^1 has the same 
distribution as = gk,2Z, where Z follows a standard unit-variance exponential distribution. Finally, 
the ergodic capacity = C'fc(7fc^i) 7fc^2' -^A' ^fc^2' ' ) whole bandwidth is equal to 

Ck = 7fc^iIE [log (1 + gkAQf^Q?)Pk,iZ)] + 7^2^ [log (l + 5fc,2n^2^)] • (3) 

Capacity is achieved if we endow the input symbols Sk{n,m) with Gaussian distribution. This 
distribution is assumed from now on. Moreover, note that does not depend on the particular subcarriers 
l^k assigned to user k, but rather on the number of these subcarriers via parameters 7^^^ and 7^2- 
Therefore, choosing some specific subcarriers rather than others has no effect on the capacity. The 
subcarriers assignment scheme reduces thus to the determination of the sharing factors 7^^, 7^2- Finally, 
the multicell resource allocation problem can be defined as follows. 

K" 

Problem 1. Minimize the power spent by the three base stations Q = X^^'^^ i^fc 1 + ^fc 2^fc 2) 

c=A,B,C k=l 

W.r.t {7^ j^,7^ 2) -ffc 1) Pk2}c=A,B,c under the following constraints: 

' k=i...k'' 

1 - 

CI :VA;,Cfc(7^,i, 7,^,2, n',i,Pfc^,2) > Rk C3 : ^^7^,2 = 



3 



K" 



C2 : 7m = « C4 : Vfc, 7^, > (i = 1, 2) . 

k=l 

As a matter of fact. Problem 1 cannot be solved using convex optimization tools. Anyhow, even if we 
were able to propose a method to solve this problem (as we did in [19] in the case of 1-D networks), 
such a method would be very costly in term of computational complexity. It is therefore of interest to 
propose practical allocation algorithms that provides suboptimal solutions to Problem 1. 
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IV. Proposed Resource Allocation Algorithm 



In [20], we showed that in 1-D cellular networks, any global solution to Problem 1 has the following 
asymptotic property: The power allocated to users who modulate both protected and non protected 
subcarriers becomes negligible as the number K of users increases. One can thus suggest the suboptimal 
(w.r.t Problem 1) resource allocation algorithm given below. For a given user k in cell c, we denote by 
{xk,yk) his/her position in the Cartesian coordinate system associated with this cell (see Figure 2). In 
our algorithm, we use a continuous function d^ubopt(-) (where D stands for the radius of 

the cell as shown in Figure 2) to define a curve that separates the users of each cell c into two subsets. 
The first subset Xj contains the users who are closer to the base station than this curve. These users 
are constrained to modulate only non protected subcarriers 3. The second subset JCp contains the rest of 
users who are constrained to the protected subcarriers Tc- 



Note that {'^siibopt(-)}c=A,B,c are fixed prior to resource allocation. Relevant selection of these curves is 
postponed to Subsection V-B. It merely relies on the asymptotic analysis carried out in Subsection V-A. 

A. Resource Allocation for Interfering Users {^^j}c=A,B,C 

For users in each cell c, resource allocation parameters in the protected subset are arbitrarily 
set to zero i.e., 1^ 2 ~ 2 ~ ^- Recall the definition of Q1 = J2k&X' ^k i^ki "-^^ average power 
transmitted by base station c {c = A, B, C) in the interference subset J. For each cell c, denote by c and 
c the other two cells. For example, A = B and A = C. Define Ck{ 7^ i, Pk v Qi' Qi) ergodic 
capacity associated with user k obtained by plugging lk2~ ^k2~^ (^)- Parameters 7^ i,-P|i for 
users in {'Xfj}c=A,B,c can be obtained as the solution to the following multicell allocation problem. 



%'j = {ke{l...K']\yk<dl,^Uxk)] 



%'p = {ke{l...K'^]\yk>dl,^Uxk)]- 



Problem 2. [Multicell problem in band 3] Minimize the total transmit power 




c=A,B,C keJCj 



' ' k=l...K'' 




CI : Vc, VA: G Xj,Rk < Ckij^,!, Pk,i,QlQl) 



C2 : Vc, 7fc, 



C3:Vc,VA:G3Cf, 7^,1, ^fc, 



> . 



k&OCj 
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Remark 2. Problem 2 may not be always feasible. Indeed, since the protected subcarriers are forbidden 
to users 3Cp the multicell interference may in some cases reach excessive levels and prevent some users 
from satisfying their rate requirements. Fortunately, we will see that if curves {d^^ubopti-)}c=A,B,c 
relevantly chosen, then the latter problem is feasible, at least for a sufficiently large number of users. 



One can use an approach similar to [19], [20] to show that any global solution to the above problem 
satisfies the following property. There exist six positive numbers {/3i, Qi}c=A,B,c (where /3f is the 
Lagrange multiplier associated with constraint C2 of Problem 2) such that: 



pc 



9k,iiQi,Qi] 
Rk 



r\9kAQlQl)/3D 



C(9fc,i(Q?,Qf)/3f) ' 
where /(.) and C(.) are increasing functions defined on by 



(4) 
(5) 



E 



l+xZ 



C{x)=E[log{l + rHx)Z)], 



(6) 



/ ^(.) being the inverse on M_|_ of /(.) w.r.t the composition of functions, and where for each c = A, B,C 
and for a fixed value of Ql and Q'j, , Q^) is the unique solution to the following system of equation: 

Rk 



E 



C{gk,iiQi,Ql)Pf] 



a , 



[9kAQl,Ql)] ' r\9k,i{QlQ\m 
Cigk^iQlQDPl) 



(7) 
(8) 



Note that equation (7) is equivalent to the constraint C2: lki~^' while equation (8) is nothing else 
than the definition of the average power Ql = X]fc6JC'= ^k i^ki transmitted by base station c in subset J. 
We now prove that when Problem 2 is feasible, then the system of six equations (7)-(8) for c = A^B,C 
admits a unique solution pf-, Qf, (3^, Qf and that this solution can be obtained by a simple 

iterative algorithm. Focus on a given cell c (c = A, B , C) and consider any fixed values Qf, Qf- Denote 
by P (Qf , Qf) the rhs of (8) i.e.. 



k&Xj 



[9k,i{Ql QD] ' r\9kAQl QDpD 

C{9kAQlQ\)Pl) 
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where is defined as tiie unique solution to (7). The value I'^ (Qi,Qi) 



can be seen as the minimum 



power that should be spent by base station c on the interference subcarriers J when the interference 
produced by base stations c and c is equal to Q1 and Q1, respectively. Since (8) should be satisfied for 
c = A, c = B and c = C, the following three equations hold 



The triple {Qf, Qf , Qi) is therefore clearly a fixed point of the vector-valued function I{Qf, Qf , Qf) 



As a matter of fact, it can be shown that such a fixed point of I is unique. This claim can be proved 
using the following lemma. 

Lemma 1. Function I is such that the following properties hold. 

1) Positivity: I(Q^,Q^,Q^) > 0. 

2) Monotonicity: If > Q^',Q^ > Q^',Q^ > then I{Q^,Q^,Q^) > I{Q^' ,Q^' ,Q^'). 

3) Scalability: for all t > 1, tI{Q^,Q^,Q^) > I{tQ'^,tQ^ ,tQ^). 

The proof of Lemma 1 uses arguments which are very similar to the proof of Theorem 1 in [23]. 
Function I is then a standard interference function, using the terminology of [24]. Therefore, as stated 
in [24], such a function I admits at most one fixed point. On the other hand, the existence of a fixed 
point is ensured by the feasibility of Problem 2 and by the fact that (9) holds for any global solution. 
In other words, if Problem 2 is feasible, then function I does admit a fixed point and this fixed point is 
unique. In the latter case, the results of [24] state furthermore that a simple fixed point algorithm (such 
as Algorithm 1 given below) applied to function I converges necessarily to its unique fixed point. 

Remark 3. Note that in Algorithm 1, the only information needed by each base station c (c = A, B,C) 
about the other two cells c, c is the current value of the powers Q\, Q\ transmitted in the interference 
band 3. This value can i) either be measured by base station c at each iteration of Algorithm 1, or ii) it 
can be communicated to it by base stations c and c over a dedicated link. In the first case, no message 
passing is required, and in the second case only few information is exchanged between the base stations. 
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l^(Qf,Qf), Qf = l''iQt,Q?), Qf = /^(Qf,Qf). 



{Qt,Qf,Q?) = iiQt,Q?,Q?)- 



(9) 
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Algorithm 1 can thus be implemented in a distributed fashion. 

Of course, the feasibility of Problem 2 depends on the choice of the separating curves {(isubopt(-)}c=A,B,c- 
Section V addresses the relevant selection of these curves such that Algorithm 1 converges for a suffi- 
ciently large number of users. 

B. Resource Allocation for Protected Users {3^p}c=7l,B,c 

Since users 3Cp in each cell c are constrained to modulate only the subcarriers of subset 7c, they are 
not subject to multicell interference. Resource allocation for such users can thus be done independently 
in each cell by solving a simple single cell optimization problem which is a special case of Problem 2. 
Focus for example on cell A. One can show [19] that the resource allocation problem for users of this 
cell is convex in variables {7fc^2' ''^A^2}feGX^' where w^2 ~ ^k2^k2- solution can be obtained by 
solving the associated KKT conditions and is given by: 

Pk,2=9klr\9ka^2) (10) 

A 

Parameter is obtained by writing that constraint = holds as the unique solution to: 

V = ^ (12) 

Resource allocation parameters for users of cells B and C can be similarly obtained. The following 
procedure performs the above resource allocation for protected users. 

C. Summary: Distributed Resource Allocation Algorithm 

The proposed distributed resource allocation scheme is finally summarized by Algorithm 3. 

D. Complexity Analysis 

By referring to Algorithm 2, it is straightforward to verify that resource allocation for protected users 
can be reduced to the determination in each cell c of the value of /Sg, which is the unique solution to 
the equation Ylk^oc^ c(g^li3'') ~ ^^^^^ function x 1— 1/C(x) is convex, the latter solution can be 
numerically obtained by any of the classical zero-finding algorithms of the convex optimization literature 
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such as the gradient method [25]. Denote by A'giad the number of iterations required till the convergence 
of such a method. Each one of these iterations requires a computational complexity proportional to the 
number of terms in the Ihs of the equation. The overall computational complexity of finding /Jg is therefore 
of order 0{Ngi-adK). In the same way, one can show that each iteration of Algorithm 1 can be performed 
with a complexity of order ©(A'grad^^)- Let A'iter designate the number of iterations of Algorithm 1 needed 
till convergence (within a certain accuracy). The overall computational complexity of Algorithm 1, and 
hence of Algorithm 3 as well, is thus of the order of ©(A^iter-^Vgrad^^)- Our simulations showed that 
Algorithm 1 converges relatively quickly in most of the cases. Indeed, no more than N^^i- = 15 iterations 
were needed to reach convergence within a very reasonable accuracy in most of the practical situations. 

V. Determination of Curves {dgy3Qp^(.)} and Asymptotic Optimality of Algorithm 3 

The aim of this section is to relevantly select the separating curves d-^^^^^{.), d^bopt(-) ^^'^ ^subopt(-)- 

For that sake, we consider the case where the number K of users tends to infinity in a sense that 

will be clear later on, and we prove Theorem 1 (see Subsection V-C) which states the following. 
There exist curves {d'^ubopti-)}c=A,B,c such that the transmit power of Algorithm 3 converges as 

K — )• oo to the limit total power of an optimal solution to the joint allocation problem (Problem 1). 
Otherwise stated. Algorithm 3 is asymptotically optimal if the separating curves are well chosen. In 

order to prove this result, we first characterize the form and the total transmit power of an optimal 

solution to Problem 1 in the special case where users of each cell are aligned on parallel equispaced lines. 

Indeed, we prove that the latter solution has the following "binary" property: In each cell c, there exists 

a curve dQc,(K) that separates users modulating uniquely protected or non protected subcarriers. Here, 

Qc,{K) ^ vector of parameters that will be specified later on and which depends on the system setting 

(including the number K of users). We show that as the number K of users tends to infinity, dga,{K) 

converges, at least for certain subsequences {K), to a curve de<- that can be characterized by solving a 

certain system of equations. The same system allows to compute the limit Qt = limK^ocQj^^ ■ Next, 

we consider the case of an arbitrary geographical distribution where users are not necessarily aligned on 

parallel lines. Eventhough the aforementioned binary property no longer holds in this general case, we 

show that the transmit power of an optimal solution to Problem 1 converges to the same limit Qt as in 

the case of aligned users. This result will suggest to relevantly select the separating curves c?subopt(-) °f 

the suboptimal allocation algorithm to be equal to the asymptotic optimal curves dQc{.). Thanks to the 

latter curve selection, we prove that the proposed allocation algorithm becomes asymptotically optimal. 
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A. Asymptotic Optimal Allocation 

The characterization of the asymptotic behaviour of an optimal solution to the joint resource allocation 
problem is performed by the following three steps. 

1) Step 1: Single Cell Resource Allocation: We first consider a particular case where users of each cell 
are aligned on equispaced parallel lines. Focus for example on cell A and define /"^ parallel equispaced 
lines (/"^ < K^) which pass through cell A and which are perpendicular to the axis BC as illustrated 
in Figure 3. Next, assign each one of these lines an index i G {1 . . . I^}. In the sequel, we denote by 
Lf C {1 . . . K^} the subset composed of the users of cell A located on the line whose index is i. Assume 
that the resource allocation parameters of users of cells B and C are fixed and recall the definition of 
Cfc given by (3) as the ergodic capacity of user k. The optimal resource allocation problem for cell A 
consists in characterizing {7^^, 7^2' Pk^2}k=i...K^ allowing to satisfy the rate requirements of 

all users k G {1, . . . , K"^}. The determination of these parameters should be done such that the power 
= Ek=i liiKi + lt2Pk,2 to be spent is minimum: 



}A _ ^^"^ r>A I r>A ...uh r^A ^,A r>A r>A 

under the following constraints: 



Problem 3. Minimize = Zk=i T^^M + ^^^fc 2 ^i^h respect to {7^^, 7^3. ^fc i' Pk,2}k=i...K- 



CI : yk,Rk <Ck C4 : 7^^ > 0,7^2 > 

fc=i 

k=l k=l 

Here, constraint C6 is a "low nuisance constraint" which is introduced to limit the interference produced 
by Base Station A. In other words, the power Qf = ^ki^ki which is transmitted by base station A 
on the subcarriers of subset J should not exceed a certain nuisance level Q. The introduction of C6 is 
a technical tool revealed to be useful in solving the multicell allocation problem later on. On one hand, 
note that Problem 3 is feasible for any a > and Q > since it has at least the following trivial solution. 
The solution consists in assigning zero power P^^ = on the subcarriers of subset J (so that constraint 
C6 will be satisfied), and in performing resource allocation only using the subcarriers of subset 7 a- On 
the other hand. Problem 3 can be made convex after a slight change of variables, as a matter of fact. 
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Therefore, any global solution to this problem is characterized by the KKT conditions. The simplification 
of these conditions is not presented in this paper due to lack of space. However, it can be done in a 
very similar way as in the case of 1-D cellular networks addressed in our previous work [19] leading to 
the following result. Resource allocation parameters of any of the subsets 'Cf of users located on lines 
i = I . . . have a "binary" separation property as the users of a 1-D cell. This property is summarized 
below. Define the following decreasing function for each x E M+: 

Z 



F{x) = E 



(13) 



l + f-\x)Z_ 

and let /32 and ^ designate the Lagrange multipliers associated with constraints C2, C3 and C6 
respectively. There exists a "pivot-position" on each line i such that users k G Lf who are farther 
than this position are uniquely assigned subcarriers from the protected subset "Pa (by setting = 0). 
Moreover, such "protected users" satisfy: 

On the other hand, users k ^ Lf who are closer to the base station than the pivot-position are uniquely 
assigned interference subcarriers from subset J (by setting = 0)- Such "non protected users" satisfy: 

?iASLfl, (^f^ft) > b.,ft) ^ (15) 

The proof of the above separation property uses Conjecture 1 in [19] which can be easily validated 
numerically. Inequalities (14) and (15) suggest the definition of a curve that geographically separates 
protected from non protected users of cell A. This can be done as follows. We write the variance pk 
of the channel gain of user k as pk = p{xk,yk) where p{x,y) models the path loss. Function p{x,y) 
is assumed to have the form p{x,y) = rj{\/ x"^ + where -^x^ + is the distance separating 

{x, y) from the base station, 77 is a certain gain and s is the path-loss coefficient. We also denote by 
5(2 (x, y) = ^^^2^-* the GNR on the protected subcarriers associated with a user at position (x, y). Note that 
for any user k, g2{xk,yk) = 5fc,2- In the same way, gi{x, y, Q', Q") denotes the GINR at position (x, y) 
if the interfering base stations are transmitting with power Q' and Q" on the interference subcarriers J. 
Using the above notation, we have gi {xk,yk, Qi, Qi) = {Qi, Qi) for each user k in cell A. Note 
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that for any {x,y), g2{x,y) = gi{x,y,0,0). Finally, for each d = /32, Q', Q", ^ we define 

g2{x,y)F{g2{x,y)P2) ■ (16) 



gi{x,y,Q',Q") ^f g,{x,y, Q',Q") ^ 
We{x,y) = ^-^ F [ ^-^^^ /?! 



1 + e V 1+e 

Due to (14), we have Wg{xk,yk) < for each protected user k i.e., for users farther from the base station 
than the pivot-position. Inversely, Wg{xk,yk) > for each non protected user k i.e., for users closer to 
the base station than the pivot-position. Therefore, function (/^(x) given below defines the curve that we 
are seeking and which geographically separates protected from non protected users of cell A: 



de{x) 



M 

2D-\x\ 
V3 



ifWe[x,i^^)<0 

'V3 



if min{We (x, , We (x, ^) } > 



(17) 



the unique zero of y i— W0{x, y) otherwise . 

Note in particular that the first two conditions of (17) hold in the case where the pivot-position at line x 
is located at the upper sector border y = |x|/\/3 or the lower sector border y = {2D — \x\) /\/3. When 
these two conditions are not satisfied, the existence of the zero of the continuous function y i— > W0{x, y) 
is straightforward due to the intermediate value theorem. The uniqueness of this zero can be proved by 
arguments already developed in the proof of Lemma 1 in [19]. Finally, we obtain the following lemma. 



Lemma 2. Assume that the users of cell A are aligned on parallel equispaced lines ( as in Figure 3 ) 
and that the power transmitted by base stations B and C on the non protected subcarriers 3 is set to 
Qf and Qi respectively. The global solution {T^fii 7Af2' -^A' ^fc^2}fc=i- -^'* Problem 3 is unique and 
is as follows. There exist three unique nonnegative numbers /3i, /32, ^ such that: 



1) For each k G Lf such that yk < dg{xk), 



p6 = ki {Qf,Q?)r'r' 



7m 



Rk 



C 



2) For each k G Lf such that yk > dg{xk), 



Pk\ = 

it = 



PC,2=9k2r\9k,2h) 



lk,2 



Rk 
CigkM 



n^2 = 



<2 = 



(18) 



(19) 
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where (3i, (32 and ^ are the Lagrange multipliers associated with constraints C2, C3 and C6 respectively, 
and where 9 = (/3i, /32, <5f , , Here, dQ{.) is the function defined by (17). 

The uniqueness of the above global solution can be proved using arguments similar to those of the 
proof of Proposition 1 in [19]. Note that due to the above lemma, there is at most one user in each subset 
who is likely to modulate both protected and non protected subcarriers. If such a "pivot-user" exists, 
then it is necessarily located on the curve do{.). Therefore, there are at most pivot-users in cell A. 

2) Step 2: From Single Cell to Multicell Resource Allocation: We now consider the problem of joint 
resource allocation (Problem 1) while still assuming that users of each cell are aligned on equispaced 
parallel lines. Recall the definition of as the subset of users of cell c located on line i {i = \ . . . P). 
The following lemma implies that any optimal solution to Problem 1 has in each cell the same form as 
the solution to the single cell problem given by Lemma 2. 

Lemma 3. Assume that the positions of users of each cell c G {A,B,C} are aligned on parallel 
equispaced lines. Any global solution {7^ 1 ' ^fe 1 ' 7fc 2' ^fc 2} "=^.3.0 to Problem 1 satisfies the following. 

' ' ' ' k=l...K'' 

=1 7fc ,1 -Pfci designate the power transmitted by base station c on the reused subcarriers 3. 
There exist nine positive numbers f]!^, ^'^}c=a,b,c ^i^ch that (18), (19) hold in each cell. 

The proof of Lemma 3 is provided in Appendix A. For each cell c G {A, B, C}, denote by c and c 
the other two cells and recall the definition of function ^^(x) given by (17) for any x G [—D,D] and 
9 G M^. Lemma 3 states that when an optimal solution to Problem 1 is applied, then there exists in each 
cell c = A, B,C a. curve d0c(.), where 9'^ = (/3f , , , ^'^), that separates protected users from non 

protected users. 

3) Step 3: Asymptotic Performance of the Optimal Resource Allocation: Denote by 9^^'^^^ = 
j32^^\ Qi^^\ Qi^^\ jtqj. c = a, B,C any set of parameters chosen such that Lemma 3 holds. 
Superscript (K) is used in order to stress the dependency of the above parameters on the number of 
users K. We now characterize the behaviour of 9"^'^^^ as the number K of users tends to infinity. Once 
the behaviour of 9^'^^^ determined, the asymptotic behaviour of both the separating curves d0a.{K) (.) and 
the total transmit power of the optimal solution to Problem 1 can be fully characterized. Assume that the 
total number K of users tends to infinity in such a way that K'^ / K — 1/3 i.e., the number of users in 
each cell is asymptotically equivalent. Denote by B the total bandwidth of the system. Define as the 
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target rate of user k in nats/s i.e., = BRk where Rk is the data rate requirement of user k in nats/s/Hz. 
Since the sum ^j^. of rate requirements tends to infinity, we let the bandwidth B grow to infinity and 
we assume that K/B — >■ t where t is a positive real number. We use in the sequel the notation /^'(^) to 
designate the number of parallel equispaced lines in cell c. We also assume that I^'(^) is such that 

jCiK) lip ,0. 

In order to simplify the proof of the results, we assume without restriction that the rate requirement for 
each user k is upper-bounded by a certain constant rmax where rmax can be chosen as large as needed. 
We also assume that for each user k, > e where e > can be chosen as small as needed. 



As a matter of fact, sequences /3^'*-^\ (32^^\ Qi^^\ Q2^^\ Sf^'^^^ are upper-bounded (refer to 
Appendix E in [22] for the proof). One can thus extract convergent subsequences from the above 
sequences. With a slight abuse of notation, ^'^'(^) = f32^^\ Qi^^\ Qi^^\ designate 

from now on these convergent subsequences and their respective limits will be denoted by 9'^ = (/3p /32, 
Qv Qv ^'^)- "^^^ provide a system of equation satisfied by the accumulation points O'^ = /Jg 
Q1, Q\, ^'^). Due to Lemma 3, the power Q^'*^^^ = Yl,k=i 7fc i-^fc i transmitted by base station c on the 
non protected subcarriers J can be written as 

'^m _ sr p, T /^^'(^) n^'(^) 



where function 3~ is defined as 

f-l ( gi{x,y,Q',Q") 

<7i(x,y,Q",Q")c(^ii^f§^/3 



^ V 1+f 

7{x, y, /3, Q', Q", i) = ^ / , /, , (21) 



for each {x, y, f3, Q', Q", ^) G [-D, -D] x [e, D] x R\. While the first term in (20) represents the power 
allocated to all the users of cell c that are uniquely assigned non protected subcarrier from subset J, the 
second term in the same equation represents the power transmitted to the (at most) I'^'^^^ pivot-users in 
the same subset. Since we assume that ^ 0, one can show [22] that the latter term is negligible 

with respect to the first term and that it tends to zero as K — oo. Thus, it will be denoted in the sequel 
by o/^(l), where Oi^(l) stands for any term that converges to zero as K tends to infinity. Define for each 
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cell c = A, B,C the following measure zy'^'(^) on the Borel sets of M+ x M x M_|_ as 

1 

J, L) = — J, L) , 



k=l 

where I, J, L are intervals of R+, R and M+ respectively and where 5{rk,Xk,yk) '^'^^ Dirac measure at point 
jXk,yk) i-^-> ^{rk,Xk,yk)(^ ^ '-^^ ^) — ^ '^^ Tk ^ I ■> Xk ^ J , Uk €i L and 6(^^^^^^^y^-^{I , J, L) — otherwise. 
Note that zy='W(/, J,L) can be interpreted as the number of users of cell c whose rate requirement in 
nats/s is inside /, whose x-coordinate is inside J and whose y-coordinate is inside L, normalized by 
A'^. In other words, measure zy^'(^) characterizes both the geographical distribution of users in cell c and 
their attribution to the different rate requirements. Replacing Rk (in nats/s/Hz) by (^^^^^^ in (20), 
we obtain 



-D Jo J~D Jmsix{\x\/V3,e} ^ ^ 

(22) 

In the sequel, we assume that the following holds. 

Assumption 1. As K tends to infinity, measure zy'^'(^) converges weakly to a measure z^^. Moreover, 
is the measure product of a limit rate distribution C,^ times a limit location distribution X^. Finally, is 
absolutely continuous with respect to the Lebesgue measure on M?. 

Note that given the definition of (^'^, the value defined as 



R'=\l rdCir) (23) 
"J JO 

represents the total average rate requirement per channel use in cell c. Here, recall that t is the limit of 
K/B as K ^ oo. It is intuitive that Ql'^^^ as given by (22) converges in this case to a constant 
defined by 

Qi = r r ^ f^' ■ (24) 

Jx=-D Jy=max{\x\/V3,e} ^ ' 

Using the same approach as above and recalling that g2{x,y) = gi{x,y,0,0), one can show that the 
power Q2^^^ transmitted by base station c on the protected subcarriers Tc converges as A — )• oo to 

D 



Ql = R' j / ^ :r(x,2/,/32^0,0,0) (iA^(x,y). (25) 

J x=—D J y=dgt: (x) 
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Now recall the expression of 7^ ^ given by Lemma 3 for all users k satisfying yk < dga,(K) (xk)- Plugging 



the latter expression into constraint C2: Yl^=i 7fc 1 = of Problem 1, we obtain 

fce{i...A'=} fce{i...KC} 



g(rE, y, Q', Q", ^ = . , „, „„, . (27) 



where we defined 



for each positive x, y, (3, Q', Q" and ^. It is thus quite intuitive that equation (26) leads as A' 00 to 

"D i-de<:{x) 
~D ^3/=max{|x|/v^,e} 

^K' r _ l-a 



[ S{x,y,(31,Ql,Ql,ndX'{x,y) = a. (28) 

Similarly, we can show that constraint C3: Yl,k=i 7fc 2 ~ °f Problem 1 leads as /^T — )• 00 to 

/ / ^ g(x,y,/32^0,0,0)dA^(x,y) = — ^. (29) 

J-DJdg.(x) -3 



Remark 4. Equations (24)-(28)-(29) characterize the asymptotic behaviour of f3'^'^^\ I32^^\ Qi^^\ 
^^'(^) in the case where users of each cell are aligned on parallel equispaced lines. The generalization 

to the case of an arbitrary setting of users is not straightforward, since Lemma 3 does not necessarily hold 

c (K) c (K) c (K) 

in this general case. Nonetheless, the lemma given below states that sequences 

/32'^"^ Q^'^"^ 

^c,{K) /jQ^g f/jg same asymptotic behaviour as given by (24)-(28)-(29) even if users are not aligned on 
parallel lines. The proof of this lemma relies on the following approach. We define in each cell a set of 
parallel equispaced lines similar to the lines in Figure 3. We next consider the projection of users positions 
on these lines using two distinct projection rules. This way, we are able to exploit equations (24)-(28)- 
(29) to solve the two resulting optimization problems. If the number of the latter lines is well chosen, then 
the perturbation of the location of each user will also be small. The optimization problem can therefore 
be interpreted as a perturbed version of the initial problem. The next step is to demonstrate that this 
perturbation of the initial setting of users does not alter the accumulation points of sequences 

can be done by properly selecting the way the number of lines scales with K. 

Lemma 4. Assume that K = K'^ + K^ — ;> 00 in such a way that K/B ^ t > and K'^/K — 1/3 
for c = A, B,C. The total power Q^"* = J2c=A B c (^fc i^k 1 + 7fc 2^fc 2) '^^y optimal solution 
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to Problem 1 converges to a constant Qt- The limit Qt has the following form: 



Qt= J2 



c=A,B,C 



( 



-D -'max{|3:,-|/v^,e} 




D J dgc (x) 



g-(x,2/,/32^0,0,0) dX^{x,y) 



(30) 



where 6'^ = /Jj, Ql, Qi,S,'^) and where for each c = A,B, C, the system of equation (24)-(28)-(29) is 
satisfied in variables 9^, Ql- Here, {x,9) i— )■ dg{x) is the function defined by (17). 

Moreover, for each c = A, B,C and for any arbitrary fixed value {Qi,Qi,Qi) = {Qi, Qf, Qf), 
the system of equation (24)-(28)-(29) admits at most one solution {(31, /32, i^). 

Lemma 4 states that the Umit Qt of the total transmit power can be computed once we have found a set 
of parameters , /Sg, Qi,i'^}c=A,B,c that satisfy (24)-(28)-(29) in the three cells c = A,B, C. However, 
these twelve parameters are underdetermined by this system of nine equations. We are nonetheless capable 
of finding {/3J, /Sg, Qi,C'^}c=A,B,c such that the above lemma holds. This can be done thanks to the fact 
that Qt is the limit of the transmit power of an optimal solution to the joint resource allocation problem. 
Therefore, {(^l, P2^QI,£,'^}c=a,b,c can be chosen as any set of parameters that satisfy the system of 
equation (24)-(28)-(29) in the three cells A, B, C and for which the total power Qt as given by (30) 
is minimal. To that end, we propose Algorithm 4 which performs an exhaustive search w.r.t points 
{Qf, Qf,Qi) inside a certain search interval. In practice, the set of points {Qf,Qf-, Qi) probed by the 
above algorithm can be determined by resorting to numerical methods. 

B. Selection of Curves {dLfo„p,(-)}c=A,B,c 

We now proceed to the relevant determination of the separating curves {c^subopt(-)}c=yl,B,c associated 
with the proposed allocation algorithm (Algorithm 3). We propose to set d'^^,^{x) such that 



where {6'^)c=a,b,c is the output of Algorithm 4 and where {x, 9) i— )■ de{x) is the function defined by (17). 

Remark 5. Note that the asymptotic separating curves dgc{.) do not depend on the particular config- 
uration of the cells, but rather on an asymptotic description of the network i.e., on the average rate 
requirement and on the asymptotic distribution of users. 



\/x G [-D,D], <ubopt(3;) = de'={x) ,c 



A,B,C . 
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Remark 6. Cun'es do>:{.) can be set before the base stations are brought into operation. They can also 
be updated once in a while if or are subject to changes. However, since such changes are typically 
slow, computational complexity of determining do>:{.) is not a major issue. 

C. Asymptotic Optimality of the Proposed Algorithm 

(K) 

Denote by Qsubopt ^^^^^ transmit power of Algorithm 3 in the case where the separating curves 
{'^subopt(-)}c=A,B,c> are selected using Algorithm 4. Recall the definition of as the total transmit 

power of an optimal solution to the multicell resource allocation problem (Problem 1). The following 
theorem states that Algorithm 3 is asymptotically optimal. Its proof is provided in [22]. 

Theorem 1. Assume that the separating curves {d'^uhopti-)}c=A,B,c <^f^ such that ~ dQc(x) 

for all X G [—D,D], where {0^)c=A,B,c die output of Algorithm 4 and where {x,9) i— )■ dg(x) is the 
function defined by (17) /or any x G [—D,D] and 6 G M^. The following equality holds: 



lim Q[1\, = lim ) = Qt 



where Qt is the constant defined by Lemma 4. 

Note that the above theorem implies that Q^^^^^i is bounded, at least for sufficiently large K. This 
means that there exists Kq g N* such that Problem 2 is feasible for all K > Kq (refer to Remark 2). 

D. Selection of the Best Reuse Factor 

During the cellular network design process, the selection of a relevant value of a allowing to optimize 
the network performance is of crucial importance. In practice, the reuse factor should be fixed prior to 
resource allocation and it should be independent of the particular cells configuration. Recall the definition 
of qi^^ = Qf\a) as the total transmit power associated with an optimal solution to the resource 
allocation problem. We define the optimal reuse factor as the value Oopt that minimizes the asymptotic 
transmit power Qria) = limA'-5.oo Qt^\ci:) (given by Lemma 4) i.e., 

Oopt = arg min Qt ■ (31) 

ae[o,i] 

In practice, Oopt can be obtained by computing QT{a) for different values of a in a grid. Note that 
computational complexity is not an issue here (refer to Remark 6). 
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VI. Numerical Results 

In our simulations, we considered the classical "free space propagation model" with a carrier frequency 
/o = 2AGHz. Path loss in dB of user k in cell c (c = A, C) is thus given by pk{dB) = 201ogio(|cA;|) + 
100.04, where \ck\ stands for the distance between user k and base station c. The thermal noise power 
spectral density is equal to Nq = —170 dBm/Hz. Denote by S the surface of any of the considered 
sectors of cells A, B, C. Each one of these three sectors is assumed to have the same uniform asymptotic 
distribution A of users, where d\{x, y) = dxdy/S. The average rate requirement R'^ in bits/s/Hz (defined 
in nats/sec/Hz by (23)) is assumed to be the same in each cell: = = R<=' = R. 
Selection of the reuse factor 

In Figure 4, we plot Oopt defined by (31) for different values of the average rate R. As expected, Oopt is 
decreasing with respect to R. Indeed, the larger the value R, the higher the level of interference, and the 
greater the number of users that should be assigned protected subcarriers. 
Selection of separating curves {d^ubopt(-) }c=a,b,c 

Once the reuse factor is set to the value aopt, the separating curves {c?subopt(-)}c=A,B,c associated with 
the proposed suboptimal allocation algorithm should be chosen to be equal to the asymptotic optimal 
curves {de<={.)}c=A,B,c given by Lemma 4. Since we are considering the case where the asymptotic 
distribution of users is the same in the three sectors. Algorithm 4 yielded in all our simulations three 
identical separating curves i.e., for all x € [—D,D], d^^^p^{x) = li^bopt(^) = '^subopt(^")- Figure 5 plots 
^subopt(-) f^'" different values of the average rate R. 
Performance of the proposed allocation algorithm 

From now on, the positions of users in each sector are assumed to be uniformly distributed random 
variables. We also assume that all users have the same target rate, and that = = . In 
the sequel, tt = 'Ylik=i '^fc designates the sum rate per sector measured in bits/s. Let us study the 
performance of the proposed allocation algorithm (Algorithm 3) in the case where the separating curves 
{'^subopt(-)}c=A,_B,c are selected as in Subsection V-B (see Figure 5). 

We first validate the asymptotic optimality of Algorithm 3. To that end, we consider 5 values of 
the number K of users comprised between 30 and 300. For each one of these values, the system 
bandwidth B = B{K) is chosen such that K/B = t = 15 x 10~^. For example, the bandwidth is 
equal to 5 MHz when K = 75 i.e., when = = = 25. This way, the number of users 
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increases in accordance with the description of the asymptotic regime given earlier in Section V-A. 

(K) 

Next, we compute the transmit powers Qsubopt spent when Algorithm 3 is applied for a large number of 
realizations of the random positions of users. We finally evaluate the associated mean value E Q[^^^l,p^ 
(expectation is taken w.r.t the random positions of users) and compare it with the asymptotic optimal 

(K) 

transmit power Qx = liuipc^oc Qt given by Lemma 4. The results of this comparison are illustrated 
in Figure 6. Note that the difference between Q^^^J^p^ and Qt decreases with the number of users. This 
difference can be considered negligible even for a moderate number of users equal to 50 per sector. This 
sustains that the proposed allocation algorithm is asymptotically optimal. 

From now on, the system bandwidth B is equal to 5 MHz and the number of users per sector is fixed to 
25. In Figure 7, we compare the proposed algorithm with the allocation scheme introduced in [23]. In the 
latter work, the authors set the value of the reuse factor a to one i.e., all the available subcarriers are reused 
in all the cells. The resource allocation problem they address consists (as in our paper) in minimizing 
the total power that should be spent by the network in order to achieve all users' rate requirements. In 
this context, they propose a distributed iterative allocation algorithm similar to Algorithm 1. The main 
difference is that, while Algorithm 1 is only applied to a subset OCf U %f U OCf of users, the algorithm 
of [23] is applied to all the users in the network. As a matter of fact, this difference has no significant 
effect on the computational complexity of the scheme of [23], which is also of order O(A^iter-^grad^) as 
Algorithm 1 (see Subsection IV-D). However, Figure 7 shows that for all the different values of the sum 
rate rx, considerable gains can be achieved by applying our resource allocation scheme instead of that 
of [23] without any additional computational complexity. 

Convergence rate of Algorithm 1 

We plot in Figure 8 the number A'^iter of iterations of Algorithm 1 as a function of the required accuracy 
i.e., the maximum relative change in the transmit powers Q1 (c = A,B,C) from iteration to another 
beyond which convergence of the algorithm is achieved. Figure 8 shows that Algorithm 1 converges 
quickly within a very good accuracy even for a sum rate as high as 9 Mbps. 

Performance of the proposed allocation algorithm in the discrete case 

We now address the so-called discrete case where the sharing factors should be integer multiples of 

In this context, we propose the following approach to compute the resource allocation parameters. We first 

apply Algorithm 3 to obtain the continuous-valued sharing factors 7^ ^ and 7^ 2 for c = A, B, C. Next, 
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we round the number of assigned subcarriers 7^ j^A^ and 7^ ^ the nearest smaller integer. In order to 
compensate for the slight decrease of each sharing factor due to rounding, the power allocated to each 
user should be slightly increased so as to keep the same achievable rate. To that end, the power allocation 
should be recomputed (this time, keeping fixed sharing factors). This can be achieved by straightforward 
adaptation of Algorithms 1 and 2. In Figure 9, we plot the required transmit power in the discrete case 
assuming the following values of the total number N of subcarriers: N = 72, 192, 360 as recommended 
in WiMax [26]. Figure 9 shows that our allocation algorithm continues to do relatively well even after 
rounding the sharing factors, provided that the total number of subcarriers is moderately large enough. 
Performance of the proposed allocation algorithm in larger networks 

We now turn our attention to the 21 -sector network of Figure 10. This network is composed of 7 duplicates 
of the 3-sector system of Figure 1. Note from Figure 10 that the subcarriers of subsets Ta, '^b, '^c are 
no more interference-free. However, the number of interferers for users modulating in these subsets is 
always smaller than the number of interferers for users modulating in subset 3. 

In this context, we propose the following procedure. We first fix the separating curve {dsubopt(-)} 
each sector c and the reuse factor a to the values given by Sections V-B and V-D respectively i.e., as 
if the network were composed of only three sectors. Note that Algorithm 2 cannot be applied anymore 
since the users outside the curves {dsubopt(-)} subject to multicell interference. Instead, we apply 

a straightforward adaptation of Algorithm 1 to the case of more than three sectors. In Figure 11, we 
plot both the transmit power of the above proposed algorithm and that of the distributed and iterative 
allocation scheme of [23] (both averaged w.r.t the random positions of users) assuming a 21 -sector 
setting. We note from Figure 1 1 that, while the scheme of [23] fails to converge for sum rates r-r larger 
or equal to 9 Mbps, our allocation algorithm converges in all the considered cases. It furthermore results 
in considerably smaller transmit powers. However, comparing Figures 7 and 1 1 reveals that the transmit 
power of the proposed algorithm is significantly larger in the 21 -sector setting than in the 3-sector setting. 
Reducing this gap requires a large amount of research and is out of the scope of this paper. 

VII. Conclusions 

In this paper, we addressed the problem of resource allocation for the downlink of a sectorized OFDMA 
network assuming fractional frequency reuse and statistical CSI. In this context, we proposed a practical 
resource allocation algorithm that can be implemented in a distributed manner. The proposed algorithm 
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divides users of each cell into two groups which are geographically separated by a fixed curve: Users of 
the first group are constrained to interference-free subcarriers, while users of the second are constrained 
to subcarriers subject to interference. If the aforementioned separating curves are relevantly chosen, then 
the transmit power of this simple algorithm tends, as the number of users grows to infinity, to the same 
limit as the minimal power required to satisfy all users' rate requirements. Therefore, the simple scheme 
consisting in separating users beforehand into protected and non protected users is asymptotically optimal. 
This scheme is frequently used in cellular systems, but it has never been proved optimal in any sense to 
the best of our knowledge. Finally, we proposed a method to select a relevant value of the reuse factor. 
The determination of this factor is of great importance for the dimensioning of wireless networks. 

Appendix A 
Proof of Lemma 3 

Notations. In the sequel, :x.abc represents a vector of multicell allocation parameters such that ^abc = 
[^a'',^b^,^bT where = [(P^)^, (7^)^]^, = [(P^)^, (7^)^]^ and xc = [(P^)^, (7^)^]^, 
and where for each c = A,B,C, P"^ = ^'i,2> • • • > ^^^2]^ and 7 = [71,1, 7i,2> • • - ih,!^ 

7^e2]^- We respectively denote by Qi(xc) = J2k^k,i^k,i Q2(xc) = Sfc 7^,2-^^,2 ^he powers 
transmitted by base station c in the interference subset 3 and in the protected subset !Pc- When resource 
allocation :x.abc is used, the total power transmitted by the network is equal to Q{'x.abc) = X]c Qi{'^c) + 
(52(xc)- Recall that Problem 1 is nonconvex. It cannot be solved using classical convex optimization 
methods. Denote by x^^^ = [x^ ,x^ ,x^ ] any global solution to Problem 1. 

Characterizing x^^^^ via single cell results. 

From ^*A^Q we construct a new vector :x.abc which is as well a global solution and which admits a 
"binary" form: for each cell c, 7^ = if > dg^{xf^) and 7^ , = if < dQ^{xk), for a certain curve 
dg^{.). For cell A, vector x^ is defined as a global solution to the single cell Problem 3 when 

a) the admissible nuisance constraint Q is set to Q = (5i(x^), 

b) the gain-to-interference-plus-noise-ratio in subset 3 is set to g^^i = gk,i (Qi(x^)) Qi(x^)). 
Vectors x^ and x^ are defined similarly, by replacing Ahy B ot C in the above definition. Denote by 
XABC = [xA^5 Xfi"'", xc^]^ the allocation obtained by the above procedure. The following claim holds. 

Claim 1. Resource allocation parameters x^bc '^f^d '^^bc coincide: li-ABC = '^^bc- 
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Proof: It is straightforward to show that 'x.abc is a feasible point for the joint multicell problem 
(Problem 1) in the sense that constraints C1-C4 of Problem 1 are met. This is the consequence of the 
low nuisance constraint (5i(xc) < (5i(x*) which ensures that the interference which is produced by 
each base station when using the new allocation 'x.abc is no bigger than the interference produced when 
the initial allocation 'x*abc i^ used. Second, it is straightforward to show that 'x.abc is a global solution 
to the multicell problem (Problem 1). Indeed, the power (5i(xc) + Q2{^c) spent by base station c is 
necessarily less than the initial power Qi(x*) + (52(x*) by definition of the minimization Problem 3. 
Thus Q{'x.ABc) < Q{'^*ABc)- course, as x^^^ has been chosen itself as a global minimum of Q, 
the latter inequality should hold with equality: Q{'x.abc) = Q{^*abc)- Therefore, x^^^ and 'x.abc are 
both global solutions to the multicell problem (Problem 1). As an immediate consequence, inequality 
(5i(xc) + Q2(xc) < <5i(x*) + Q2(x*) holds with equality in all the three cells c = A,B, C: 

Ql(Xc) + Q2(Xe) = Qli^l) + Q2{0 . (32) 

Clearly, x^ is a feasible point for Problem 3 when setting Q = Qi(x^) and g^ i = (?a;,i(Qi(x^), 
(5i(x^)). Indeed constraint C6 is equivalent to Qi{x.\) < Q and is trivially met (with equality) by 
definition of Q. Since the objective function Qi(x^) + Q2{'x*a) coincides with the global minimum as 
indicated by (32), x^ is a global minimum for the single cell Problem 3. By Lemma 2, this problem admits 
a unique global minimum :ka- Therefore, x^ = x^i. By similar arguments, x^ = x^ and x^ = x^. ■ 
We thus conclude that any global solution ^*abc multicell Problem 1 satisfies equations (18), 

(19), where gf^^i in the latter equations coincide with gj^ i = gj^ i ((5i(xg), (5i(xS)), and where for each 
cell c G {A, B, C}, c and c denote the other two cells. The proof of Lemma 3 is thus complete. 
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Figure 1. 3-cells system model and the frequency reuse scheme 




Figure 2. Fixed separating curve in cell A 
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Algorithm 1 Ping-pong algorithm for three interfering cells 

Initialization: ^ 0, Qf ^ 0, Qf ^ 
repeat 

il^t^Qi) ^ Solve (7)-(8) for c = A 

(/3f ,Qf ) ^ Solve (7)-(8) for c = B 

{Pi,Qi) ^ Solve (7)-(8) for c = C 
until convergence 
for all c = A,B,CAo 

{7,V,n?i}fcexj ^(4)-(5) 
end for 

return {7^,1, i=fc,i}c=A,B,c, fceocj 



Algorithm 2 Resource allocation for protected users 
for all c = A,B,C Ao 
/3f ^ Solve (12) 
for all ke%%Ao 

Pk,2 ^ (10) 

7^,2^(11) 
end for 

end for 

return {7fc,2' ^fc,2}c=A,i3, feexj. 



Algorithm 3 Proposed resource allocation algorithm 
for allc = ^,B,Cdo 

X'i^{kG{l...K-}\yk<di^,^^^ixk)} 
end for 

{7fc,i,-Pfc,i}c=A,B,c,fcexj ^ Algorithm 1 
{7fc,2'-Pfc,2}c=A,B,c,fceX|, ^ Algorithm 2 
return {lli,Pk,i,lk,2^Pk,2}c=A,B,c,k=i...K^ 




Figure 3. Definition of subsets {Cf}^^^ ja and of the curve dg{.) 
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Algorithm 4 Determination of {9'^}c=A,B,c 

for all iQf,Qf,Qf) do 
for c = A, B,C do 

if (24)-(28)-(29) admits a solution then 

(/3J,/3^,C'=) ^ unique solution to (24)-(28)-(29) 

2D-\x\ 

Q'^Ql + /f^ Xi^^^) J(rE, y, /3f , 0, 0, 0) d\%x, y) 
else 

Q"^ ^ oo 
end if 
end for 

QT{QtQtQi)^T.c=A,B,cQ' 

end for 

(Qf , Qf , Qf ) ^ argmin^Q^ (3r(0i^, Qf , Of) 

forc = A,S,Cdo 

end for 

return 9^,9^,9^' 



D=500 m 




R (bits/s/Hz) 

Figure 4. Optimal reuse factor vs. average rate of a sector 
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R= 1 bit/ s/ Hz 
R=IA bit/ s/ Hz 
^=1.8 bit/ s/ Hz 
7^ = 2.4 bit/ s/ Hz 



Figure 5. Optimal separating curve dgA{.) 



D=500 m, r^=7 Mbps 
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Figure 6. E ^Qfilbop, — Qt^ /Qt vs. number of users per sector 
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Figure 7. Comparison between the proposed allocation algorithm and the scheme of [23] for = = K'^ = 25 
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Figure 8. Number of iterations of Algorithm 1 vs. relative accuracy 
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Figure 9. Transmit power of the proposed algorithm in case the sharing factors are integer multiples of 
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Figure 11. Comparison between the proposed allocation algorithm and the scheme of [23] in the case of 21 sectors for K'^ = 25 



September 14, 2010 



DRAFT 



